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Abstract 


The Ebola virus is currently one of the most virulent pathogens for 
humans. The latest major outbreak occurred in Guinea, Sierra Leone 
and Liberia in 2014. With the aim of understanding the spread of in¬ 
fection in the affected countries, it is crucial to modelize the virus and 
simulate it. In this paper, we begin by studying a simple mathematical 
model that describes the 2014 Ebola outbreak in Liberia. Then, we use 
numerical simulations and available data provided by the World Health 
Organization to validate the obtained mathematical model. Moreover, 
we develop a new mathematical model including vaccination of individ¬ 
uals. We discuss different cases of vaccination in order to predict the 
effect of vaccination on the infected individuals over time. Finally, we 
apply optimal control to study the impact of vaccination on the spread 
of the Ebola virus. The optimal control problem is solved numerically 
by using a direct multiple shooting method. 

1 Introduction 

Ebola is a lethal virus for humans. It was previously confined to Central 
Africa but recently was also identified in West Africa [X]. As of October 8, 
2014, the World Health Organization (WHO) reported 4656 cases of Ebola 
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842792, 9 pp. See http://dx.doi.org/10.1155/2015/842792 Submitted 26-Dec-2014; 
revised 27-Feb-2015; accepted 28-Feb-2015. 
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virus deaths, with most cases occurring in Liberia |2j- The extremely rapid 
increase of the disease and the high mortality rate make this virus a major 
problem for public health. Typically, patients present fever, malaise, ab¬ 
dominal, headache, and asthenia. One week after the onset of symptoms, a 
rash often appears followed by haemorrhagic complications, leading to death 
after an average of 10 days in 50%~90% of infections [3H6j. Ebola is trans¬ 
mitted through direct contact with blood, bodily secretions and tissues of 
infected ill or dead humans and non-human primates um- 

The main aim of this work is to understand the dynamics of the Liberian 
population infected by Ebola virus in 2014, by using an appropriate mathe¬ 
matical model. More precisely, we begin by considering a system of ordinary 
differential equations to describe the 2014 Ebola outbreak, which is nothing 
else than an epidemic SIR model, that is, a model based on the division of 
the population into three groups: the Susceptible, the Infected, and the Re¬ 
covered mm- We simulate the model obtained by parameters estimated 
on November 4, 2014, by Kaurov m- Our objectives are to better under¬ 
stand the outbreak and to predict the effect of vaccination on the infected 
individuals over time. 

In order to deal with the epidemic of Ebola, governments have decided 
to implement tough measures. For example, some have applied quarantine 
procedures while others have opted for mass vaccination plans as a precau¬ 
tion to the epidemic mm- In this context, we consider an optimal control 
problem to study the effect of vaccination on the spread of virus. 

The text is organized as follows. In Section [2] we introduce a basic math¬ 
ematical model to describe the dynamics of the Ebola virus during the 2014 
outbreak in Liberia. After the mathematical modelling, we use the obtained 
model to simulate it in Section [3] with the parameters estimated from recent 
statistical data based on the WHO report of the 2014 Ebola outbreak [18 j. 
Following [IB! > i n Section [I] we study different cases of vaccination of indi¬ 
viduals by adding a vaccination term to the model. Section [5] presents an 
optimal control problem, which we use to study the impact of a vaccination 
campaign on the spread of the virus. We end with Section [6] of conclu¬ 
sions and future work, where the results are summarized and some research 
perspectives presented. 

2 The Basic Mathematical Model 

The objective of this section is to describe mathematically the dynamics 
of the population infected by the Ebola virus. The dynamics is described 
by a system of differential equations. This system is based on the common 
SIR (Susceptible-Infectious-Recovery) epidemic model, where the popula¬ 
tion is divided into three groups: the susceptible group, denoted by S(t), 
the infected group, denoted by I(t), and the recovered group, denoted by 
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Figure 1: The scheme of the Susceptible-Infected-Recovered model 


R(t). The total population, assumed constant during the short period of 
time under study, is given by N = S(t) + I(t) + R(t). 

To create the equations that describe the population of each group along 
time, let us start with the fact that the population of the susceptible group 
will be reduced as the infected come into contact with them with a rate of 
infection /3. This means that the change in the population of susceptible is 
equal to the negative product of (3 with S (t) and I (t): 

dS(t ) S(t)I(t) 

dt 'N ' 1 ’ 

Now, let us create the equation that describes the infected group over time, 
knowing that the population of this group changes in two ways: 


(i) people leave the susceptible group and join the infected group, thus 

adding to the total population of infected a term j3S{t)I[t)\ 

(ii) people leave the infected group and join the recovered group, reducing 

the infected population by 


This is written as 


dl(t) 

dt 


3 smt ) 

P N 


nl(t). 


( 2 ) 


Finally, the equation that describes the recovery population is based on the 
individuals recovered from the virus at rate /r. This means that the recovery 
group is increased by [i multiplied by I(t): 


dR(t) 

dt 


til(t). 


( 3 ) 


Figured] shows the relationship between the three variables of our SIR model. 
As shown in the next section, this simple model describes well the 2014 Ebola 
outbreak in Liberia for suitable chosen parameters /3 and n and appropriate 
initial conditions 5(0), 1(0) and R( 0). 
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3 Numerical Simulations of the Basic Model and 
Discussion 


With the aim to approximate the 2014 Ebola outbreak in Liberia, we now 
simulate the complete set of equations JT]) (0 that describe the SIR model: 


' dS(t) 
dt 

dm 

dt 

dR(t ) 
' dt 


p S(t)m 

N 

p S(t)m 

N 

= 


ni(t), 


(4) 


where (3 is the rate of infection and the rate of recovery. The data for the 
outbreak in Liberia is based on the cumulative numbers of total reported 
cases (confirmed and probable) taken from the official data reported by 
the Ministry of Health of Liberia up to October 26, 2014, as provided by 
the World Health Organisation (WHO) [18] (in particular, the Ebola virus 
disease cases reported each week from Liberia are given in Figure 2 of USD- 
The model fits well the reported data of confirmed (means infected) cases in 
Liberia. Note that the exact population size N does not need to be known 
to estimate the model parameters as long as the number of cases is small 
compared to the total population size [20]. To estimate the parameters of the 
model, we adapted the initialisation of S and I with the reported data of m 
by fitting the actual data of confirmed cases in Liberia. The data that we 
are studying is the number of individuals with confirmed infection between 
July 14 and October 12, 2014. The result of fitting is shown in Figure [2] The 
comparison between the curve of infected obtained by our simulation and 
the reported data of confirmed cases by WHO shows that the mathematical 
model ([ID fits well the real data by using (3 = 0.000318 as the rate of infection 
and /i = 0.0175 as the rate of recovery. The initial numbers of susceptible, 
infected, and recovered groups, are given by 5(0) = 460, 1(0) = 12, and 
R(0) = 0, respectively. The choice of 5(0) is in agreement with the data 
shown in Figure 2 of m, as the biggest value in the histogram is 460. To 
interpret the simulation results, we describe the evolution of each group of 
individuals over time and the link between them. Figure [3] shows that the 
susceptible group begins to plummet and, simultaneously, the number of 
infected begins to rise. This is due to how infectious the Ebola virus is. 

When we compare our results with those of Astacio et al. [21], when 
they modeled the 1976 outbreak of Ebola in Yambuku, a small village in 
Mongala District in northern Democratic Republic of Congo (previously 
Zaire), and the 1995 outbreak of Ebola in Kikwit, in the southwestern part of 
the Democratic Republic of Congo, we note that the evolution of each group 
obtained in our simulations is reasonable. The main difference between the 
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Figure 2: Graph of Infected I(t) obtained from ((H) with j3 = 0.000318, 
/r = 0.0175, 5(0) = 460, 1(0) = 12, and 5(0) = 0 versus the real data of 
confirmed cases for the 2014 Ebola outbreak occurred in Liberia from July 
14 (day t = 0) and October 12, 2014 (day t = 90). 



Figure 3: The evolution of susceptible (5(f)), infected (/(f)), and recovered 
( R(t )) for the period July 14 (t = 0) and October 12, 2014 (t = 90), given 
by the Ebola Liberia model (HD with f3 = 0.000318, /i = 0.0175, 5(0) = 460, 
1(0) = 12, and -5(0) = 0. The real data of confirmed cases for the 2014 
Ebola outbreak occurred in Liberia, taken from |18j , is shown with circles. 
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2014 outbreak of Ebola in Liberia and the previous outbreaks of Ebola in 
Yambuku and Kikwit studied in [21] is at the end of the evolution of infected: 
while in the previous 1976 and 1995 outbreaks of Ebola in Yambuku and 
Kikwit the infected group goes to zero, now in the 2014 outbreak of Ebola in 
Liberia the infected group remains important with a high number of infected 
individuals by the virus, which describes the reality of what is currently 
happening in Liberia. One can conclude that the 2014 Ebola outbreak of 
Liberia is much more dangerous than the previous Ebola outbreaks. 


4 Modelling and Simulation of Ebola Outbreak 
with Vaccination 

In this section, we simulate the model by using parameters estimated on 
November 4, 2014, by Kaurov m, who studied statistically the recent data 
of the World Health Organisation (WHO) [18]. In his statistical study, 
Kaurov studied the outbreak by modeling it with Wolfram’s Mathematica 
language. Here we use his SIR model for modelling and time-discrete equa¬ 
tions for resolution of the set of equations through Matlab. 


4.1 The Basic Model with Kaurov’s Estimation of Parame¬ 
ters 

Let us recall the complete set of equations considered by P3j : 


' dS(t) 
dt 

dm 

dt 

dR(t ) 
- dt 


= - l3S(t)I(t ), 

= mm - 

= vi a), 


(5) 


where, as before, /3 is the rate of infection and /i is the rate of recovery. 
The parameters obtained by Kaurov, knowing that 95% of population is 
susceptible and 5% of population is infected, are (3 = 0.2 and /r = 0.1 [13]. 
In agreement, the initial susceptible, infected, and recovered populations are 
given, respectively, by 5(0) = 0.95, 1(0) = 0.05, and R(0) = 0 (note that in 
contrast with Sections [2] and [3] the variables are now dimensionless). 

Our numerical simulations were done with the ODE solver of Matlab. 
The results are shown in Figure [4] To interpret the results, we describe the 
evolution of each group over time and the links between them. Figure 3] 
shows that the susceptible group begins to plummet due to how infectious 
the virus is and, simultaneously, the infected group’s number begins to rise. 
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Figure 4: Solution to the 2014 outbreak Ebola model © with (3 = 0.2, 
H = 0.1, 5(0) = 0.95, 7(0) = 0.05, and R( 0) = 0. 


An important dimensionless quantity that plays a key role in any SIR 
model is the well-known basic reproduction number Rq [22] . This number 
tells us how fast the disease will spread at the epidemic. For model j5|), 
the basic reproduction number is simply given by /3//i. Let us now discuss 
the comparison between the different values of Rq obtained recently. On 
September 2, 2014, Althaus studied the case of the 2014 Ebola outbreak oc¬ 
curred in Guinea, in Sierra Leone, and in Liberia [2D] • In his study, Althaus 
used the WHO data occurred in the three countries. He found different 
values for the basic reproduction number, depending on the country: in 
Guinea T?o = 1-51, in Sierra Leone the value of Rq is given by 2.53, while in 
Liberia Rq is equal to 1.59. Let us now compare between these values and 
the value for the basic reproduction number studied by Kaurov, knowing 
that his study, which was done on November 4, 2014, is more recent than 
the one of Althaus, done September 2, 2014. As the parameters estimated 
statistically by Kaurov are for the three countries, our comparison will be 
between Ro,a, the mean of the three basic reproduction numbers found by 
Althaus, and Ro,k, the number found more recently by Kaurov. The value 
of 7 ?o,a is given by 1.876, where the more recent Rq,k is given by 2. We 
claim that the basic reproduction number increased between September and 
November 2014. 
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Figure 5: The scheme of the Susceptible-Infected-Recovered model with 
vaccination ©). 

4.2 Dynamics in Case of Vaccination 

Infectious diseases have tremendous influence on human life. In last decades, 
controlling infectious diseases has been an increasingly complex issue. A 
strategy to control infectious diseases is through vaccination. Now our idea 
is to study the effect of vaccination in practical Ebola situations. The case of 
the French nurse cured of Ebola with the help of an experimental vacine is a 
proof of the possibility of treatment by vaccinating infected individuals |23| . 

In this section, we improve the SIR model described in Section T4. 11 The 
system of equations that describe the Susceptible-Infectious-Recovery (SIR) 
model with vaccination is: 

= - PS(t)I(t ) - vS(t), 

< ^ = ( 6 ) 

—= fiI(t) + vS(t), 
at 

where /3 is the rate of infection, /r is the rate of recovery, and v is the 
percentage of individuals vaccinated every day [23] • Figure 0 shows the 
relationship between the variables of the SIR model with vaccination (|6|) . 


4.3 Numerical Simulation in Case of Vaccination 

We simulate the SIR model with vaccination © in order to predict the 
evolution of every group of individuals in case of vaccination. Our study 
is based in the test of different rates of vaccinations and their effect on 
the curve of every group. The results are shown in Figures [6] and [TJ 
These figures present different rates of vaccination and their effect on the 
peak of the curve of infected individuals along time. The time-dependent 
curve of infected individuals shows that when the percentage of vaccinated 
individuals increases, the peak of the curve of infected individuals is less 
important and the period of infection (the corresponding number of days) 
is shorter while the number of recovered gets more important. In fact, in 
case of infection without vaccination, the curve of the infected group goes 
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(a) v = 0 



(b) v = 0.005 



(c) v = 0.01 

Figure 6: A solution to the SIR model © with j3 = 0.2, ^ = 0.1, 5(0) = 
0.95, 7(0) = 0.05, 7?(0) = 0, and different rates of vaccination: v = 0 (a), 
v = 0.005 (b), and v = 0.01 (c). 
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(a) v = 0.02 



(b) v = 0.03 



(c) v = 0.06 

Figure 7: A solution to the SIR model © with (5 = 0.2, // = 0.1, 5(0) = 0.95, 
1(0) = 0.05, i?(0) = 0, and different rates of vaccination: v = 0.02 (a), 
v = 0.03 (b), and v = 0.06 (c). 
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Rate of vaccination 

Maximum percentage 
of infected 

Days to the peak of 

m 

0 

17% 

26.52 

0.005 

15% 

23.47 

0.01 

13% 

22.07 

0.02 

10% 

17.09 

0.03 

9% 

13.68 

0.06 

7% 

9.28 


Table 1: Maximum percentage of infected individuals 7 corresponding to 
different rates v of vaccination and respective number of days for the peak 
of 7. 


to zero in about 90 days (see Figure [6]), where in the presence of vaccination 
the curve of the infected group goes to zero in about 45 days (see Figure [7]). 
This shows the efficiency of using vaccination in controlling Ebola. Table [T] 
presents the maximum percentage of infected individuals, corresponding to 
different rates of vaccinations, and the respective number of days for the 
peak of infected. 

5 The Optimal Control Problem 

Optimal control has been recently used with success in a series of epidemi¬ 
ological problems: see Q3E5H27] and references therein. In this section, 
we present an optimal control problem by introducing into the model (JSJ) a 
control u(t) representing the vaccination rate at time t. The control u(t) is 
the fraction of susceptible individuals being vaccinated per unit of time. We 
assume that all susceptible individuals that are vaccinated are transferred 
directly to the removed class. Then, the mathematical model with control 
is given by the following system of nonlinear differential equations: 

C -^p- = -PS(t)I{t ) - u(t)S(t), 

< ^ = / 3 ( 7 ) 

= fiI(t) + u(t)S(t), 

where 5(0) > 0, 7(0) > 0, and 72(0) > 0 are given. Our goal is to reduce 
the infected individuals and the cost of vaccination. Precisely, our optimal 
control problem consists of minimizing the objective functional 

J (u) = J I{t) + ^-u 2 (t)dt, ( 8 ) 
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Figure 8: The function I(t) of infected individuals with and without control. 


subject to 0 < u(t) < 0.9 for t £ [0, t/], where the parameters A > 0 
and tf denote, respectively, the weight on cost and the duration of the 
vaccination program. We present the numerical solution to this optimal 
control problem, comparing the results with the numerical simulation results 
obtained in Section H] via Matlab. As in Section U the rate of infection is 
given by /3 = 0.2, the recovered rate is given by n = 0.1, and we use 
5(0) = 0.95, 1(0) = 0.05, and R( 0) = 0 for the initial number of susceptible, 
infected, and recovered populations, respectively (i.e., at the beginning 95% 
of population is susceptible and 5% is already infected with Ebola). Note 
that when u(t) = 0 model ( 0 ) reduces to ( 0 ); and when u(t) = v model ( 0 ) 
reduces to ©. 

The numerical simulation of the system without control (u(t) = 0) was 
done with Matlab in Section [4.11 For our numerical solution of the optimal 
control problem, we have used the ACADO solver [28], which is based on 
a multiple shooting method, including automatic differentiation and based 
ultimately on the semidirect multiple shooting algorithm of Bock and Pitt 
[29J. ACADO is a self-contained public domain software environment written 
in C++ for automatic control and dynamic optimization. Figures [S] and 0 
show the significant difference in the number of infected individuals with and 
without control. In fact, in case of optimal control, the function I decreases 
greatly: the maximum value on the infected curve I under optimal control 
is 5.6% against 17.9% without any control (see Figure [9]). 

In Figure [101 we see that the number of susceptible 5 decreases more 
rapidly during the vaccination campaign. It reaches 11% at the end of the 
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X: 26.27 
Y: 0.179 




Figure 9: Comparison between the peaks of infected individuals I(t) in case 
of optimal control versus without control. 
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Figure 10: Susceptible individuals S(t) in case of optimal control 
without control. 



Figure 11: Recovered individuals R(t) in case of optimal control 
without control. 


versus 


versus 


14 
















0.25 


Control 


0.2 - 


0.15 


0.1 - 


0.05 


0 10 20 30 40 50 60 70 80 90 100 

Time [days] 


Figure 12: The optimal control u for problem (j7]) -(|5 |) . 


campaign against 19% in the absence of optimal control. 

Figure [TT] shows that the number of recovered individuals increases rapidly. 
The number R(tf ) at the end of the optimal control vaccination period is 
88.7% instead of 80.5% without control. Figure [12] gives a representation of 
the optimal control u(t). 

6 Conclusion 

We investigated a mathematical model that provides a good description of 
the 2014 Ebola outbreak in Liberia: the model fits well the data of con¬ 
firmed cases in Liberia provided by the World Health Organization. The 
evolution of infected individuals shows that the infected population remains 
important, describing the actual evolution of the virus in Liberia. This evo¬ 
lution is different than Ebola’s outbreak in the last decades. Our study of 
the SIR model with vaccination shows that vaccination is a very efficient 
factor in reducing the number of infected individuals in a short period of 
time and increasing the number of recovered individuals. We also discussed 
an optimal control problem related to the impact of a vaccination strategy 
on the spread of the 2014 outbreak of Ebola in West Africa. The numer¬ 
ical simulation of both systems, with and without control, shows that an 
optimal control strategy greatly helps to reduce the number of infected and 
susceptible individuals and increases the number of recovered individuals. 

As future work, we plan to include in our study other factors. In par¬ 
ticular, we intend to include in the mathematical model the treatment of 
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infected individuals with a quarantine procedure. Another interesting line 
of research is to investigate the effect of impulsive vaccination. 
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